^ ■ ELLIPTIC FLOW IN A FINAL STATE INTERACTION MODEL 

O' 

o 



< 



> 

OO 



o 



a 



A. CAPELLA("), E. G. FERREIRO ^'') 

^°' Laboratoire de Physique Theorique, Universite de Paris XI, F-91405 Orsay Cedex, France 

^^' Depto. de Fisica de Particulas, Universidade de Santiago de Compostela 

E- 15782 Santiago de Compostela, Spain 



Cn ■ We propose a final state interaction model to describe the fixed pt suppression of the yield 

of particles at all values of pr- We make an extension of the model to the motion in the 
transverse plane which introduces a dependence of the suppression on the azimuthal angle 9r. 
We obtain values of the elliptic flow V2 close to the experimental ones for all values of pr- 



"^ ! 1 Introduction 

o 



The RHIC data show azimuthal anisotropy in the production of particles in heavy ion collsions. 
This asymmetry is currently named elliptic flow V2- Two main interpretations has been applied 



o ■ in order to explain the data. In the framework of the hydrodynamical models, it is considered 

as a signal of early thermalization (Quark Gluon Plasma formation): particles tend to go in the 
Q \ direction of the strongest pressure gradients, hence preferably in the collision plan. On the other 

1-^ ' hand, at larger transverse momenta, measurements of azimuthal anisotropy are also relevant to 

the observation of jet quenching - final state interactions. In the second scenario, due to the 
k><( \ asymmetry of the overlap region of the two nuclei, the average amount of matter traversed by a 

;h ' parton depends on its azimuthal direction with respect to the reaction plane, which leads to an 



azimuthal anisotropy of the emitted jets. Our mechanism consists on a final state interactions 
approach applied to the whole pT region. In this approach, we can reproduced the large px 
suppression and the elliptic flow data. In this way, the V2 in our approach becomes medium 
dependent, while in other approaches the V2 corresponds to the gradient of the medium. 

2 The final state interaction model 

The interaction of a particle or a parton with the medium is described by the gain and loss 
differential equations which govern final state interactions: 

T -r- = ^(^ke Pk pe -^ (^ik Pi Pk , (1) 

where pi = dN ^^'^{b)/dyd?s are transverse densities and aij are the final state interaction 
cross-sections. The first term of (^ describes the gain in type i particle yield resulting from 
the interaction of k and i. The second one corresponds to the loss of type i particles resulting 
from its interaction with particle k. Consider a vr^ produced at fixed pT interacting with the 
hot medium. In the interaction, with cross-section a, the t[^ suffers a decrease in its transverse 



momentum with a pT-shid Spx- This produces a loss in the vr^ yield in a given px bin. There 
is also a gain resulting from vr'^'s produced at px + SpT and that are shifted to smaller values of 
Pt- The gain and loss diferential equation for pions is then 






-0- Pmedium PTro{b, S,y,PT) + CT Pmedium PTro{b, S,y,PT + Spr) 



(2) 



Due to the steep fall-off of the pt spectrum the loss is larger than the gain, resulting in a 
net suppression of the n^ yield at a given pT- Our equations have to be integrated between 
initial time tq and freeze-out time rj. The solution depends only on the ratio tj/tq. We use 
the inverse proportionality between proper time and densities, tj/tq = p{b,s,y)/ppp{y), where 
Ppp{y) corresponds to the density per unit rapidity for pp collisions at \/s = 200 GeV= 2.24 
fm~^ and p{h, s, y) is the density produced in the primary collisions. Our densities can be either 
hadrons or partons. In fact, at early times, densities are very high and hadrons not yet formed, 
so our equations describe final state interactions at a partonic level. At later times we have 
interactions of full fledged hadrons, and, thus, a represents an effective cross-section averaged 
over the interaction time. 

Integrating eq. (jSj) from tq to r/ and taking p{h^y^pT)-nO = dN^o/dbdydpT, we obtain the 
suppression factor 6*^0(6, y,pr) of the yield of vr'^'s at given pT and at each impact parameter, 
due to its interaction with the dense medium: 



'S'^o {y,pT,b) = exp{-a p{b, y) 
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(3) 



When 5pT tends to co, the gain term vanishes, and the survival probability has the same 
expression as in the case of J/V' suppression without cc recombination . If Spx is equal to 0, 
the loss and gain terms are identical and the survival probability becomes one. 



3 Numerical results 

In order to perform numerical calculations, we need the px distribution of the vr^'s. We use the 
following parametrization for the ratio R\j^{b,pT) in the absence of final state interactions: 
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(4) 



where PQ{b) = (n — 3)/2 < px >b, n = 9.99 and < px >b is the experimental value of < pT > at 
each b. We have also tried different parametrizations'=' for i?^^(6,pj'). Our final result depends 
little on the form of -R^^ taking a pT-shift Spr of the form p^/C. To the value R\j<^ we apply the 
correction due to the suppression factor 5^o, RAA{b-,PT) = R\A^b.,PT) S^o{b,pT)- Our results 
are shown in Fig. 1. The dashed lines correponds to the result obtained using eq. (jlj) for 
R\uAu ^'^'^ ^ pT-shift given by 5pT = p^ /(20 GeV^'^). The continous lines are obtain with 
RAA{b,PT > 5 GeV) = 1 and the pT-shift given by 6pT = Pt^ /{20 GeV^/^) for pr < 2.9 GeV 
and 6pT = p^^/{9.5 GeV^/^) ^^ ^^ > 2.9 GeV. 
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Figure 1. Left: The 7r° suppression factor RAuAu{PT,b) at ^/s — 200 GeV in the centrality bin 0-10 % 

and the corresponding RauAu(PTj b) in the absence of final state interaction. 

Right: Centrahty dependence of RauAu for pt > ^ GeV/c. The data's are from PHENIX. 



4 Elliptic flow 

Our final state interaction model takes into account the longitudinal expansion with no consid- 
eration for the motion in the transverse plane. Elliptic flow, on the contrary, results from an 
asymmetry in the azimuthal angle, so the motion in the transverse plane plays a fundamental 
role. We propose a simple extension of the model™ taking into account the different path length 
of the particles or partons in the transverse plane for each value of its azimuthal angle Or — 
measured with respect to the reaction plane. At y* ~ 0, the path length -Re^, measured from 
the center of the interaction region (overlap of the colliding nuclei) is given by 



Re^{b) = RA 
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where Ra = 1-05 A^'^ fm is the nuclear radius and sin a 
the following replacement in eq. © 
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(5) 
bsm9R/2RA- Our ansatz consists in 



(6) 



This is motivated by the fact that the duration of the interaction, as well as the density of 
the medium traversed by the particles, are expected to be proportional to the path length i?^^ 
inside the overlap region of the colliding nuclei. Due to the division by < Rg^^ > in eq. ©, the 
results of section 3 are unchanged. With the replacement @, the survival probability becomes 
^R-dependent, and the elliptic flow can be obtained as 
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Clearly, when the particle or parton moves along the reaction plane 9r = 0° its path length 
will have its minimal value and the survival probability its maximal one. On the contrary, for 
9r = 90° the path length will be maximal and the survival probability minimal. 

Using the same value a = 1.3 mb of the final state interaction cross-section and the same 
PT-shiit introduced in section 3 in order to describe the experimental values of RAA{b,PT), we 
obtain the values of V2 versus px and centrality shown in Fig. 2. 





Pt (GcV/c) 

Figure 2. Left: V2 vs. px at centrality 13 %-26 %. Right: i;2 vs. the number of participants for different 
values of pt: 0.4 (lower), 0.75 (middle) and 1.35 GeV (top). Datal^ PHENIX (black), STAR (open). 
We find a good agreement with data for the pT and the centrality dependence. For the mass 
dependence, in Fig. 3, our V2 of mesons falls below that of baryons for pT > 2 GeV, in agreement 
with data, while the hydrodynamical model predicts the same mass ordering for V2 at all px- 
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Figure 3. Left: Values of RAuAu{b, Or) as a function of the azimuthal angle On for pr > 4 GeV in various 
centrality bins. Right: V2 vs. px of differents particles for mb collisions. Data'^'are from PHENIX 

We have proposed a final state interaction model which takes into account the different 
path length of a particle in the transverse plane for each value of its azimuthal angle. In our 
approach, the mechanism responsible for the large pT suppression gives a contribution to V2- 
Elliptic flows very close to the experimental ones are obtained in the whole px region. Although 
this contribution to V2 results from an asymmetry in the azimuthal angle, it can be qualified as 
non-flow: the mechanism from which it arises (fixed px suppression) is maximal at zero impact 
parameter and thermalization is not needed. We do not claim that our mechanism gives the 
only contribution to the elliptic flow. In our opinion, our knowledge of the dynamics of the 
nuclear interaction is not sufficient to disentangle all the interpretations regarding the V2. 
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